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^ ■ Dynamical symmetry breaking is studied in an Eq GUT model of a single 

\ generation of fermions with strong 4-fermi interactions. The effective potential is 
analyzed analytically by the help of Michel's conjecture and the result is con- 
firmed numerically. We find that the Eq symmetry is spontaneously broken either 
to F4 or to 5*^(8) or G2 or S'f/(3), depending on which of the 4-fermi coupling 
constants G27 and 6*351 in the 27/351 channels is stronger. The possibilities for 
obtaining other type of breaking patterns are also discussed. 



1. Introduction 



Despite remarkable successes of the standard model based on SU (3) x SU (2) x 
f/(l), many physicists believe that there exists a more fundamental theory beyond 
it. The strongest evidences for such theories, usually called grand unified theories 
(GUT's) , are the facts that the quark/lepton charges are quantized and that 
anomaly is cancelled miraculously between quarks and leptons. In the usual sce- 
nario of GUT's, however, the spontaneous symmetry breaking required there is 
discussed by introducing some elementary Higgs fields just in the same manner as 
in the standard model. Then quite a large arbitrariness appears in, e.g., which 
representations and how many we introduce as the Higgs fields. Moreover this 
introduces too many arbitrary parameters, even more than in the standard model, 
in the Higgs Yukawa- and self-couplings. 

[3—81 

Dynamical symmetry breaking scenario is very attractive in this respect. 
There one supposes that there exist only matter fermion fields belonging to some 
representation of a gauge group G and the gauge fields of that group. Then the 
Lagrangian is uniquely determined by the gauge symmetry alone when we require 
the renormalizability (and if the fermions are all chiral). The usual Higgs fields 
arc supplied as bound states of the fundamental fermions which are formed by 
the gauge interaction dynamics itself. So, which types of Higgs fields appear is 
determined dynamically and all the parameters concerning the Higgs fields, which 
are arbitrary in the usual scenario, becomes in principle calculable. 

Even when it is difficult to solve fully the dynamics, the dynamical symme- 
try breaking scenario can give several constraints on the possible models for the 
GUT's, e.g., on possible GUT groups and/or matter contents. For instance, as 
was emphasized by Barbieri and Nanopoulos^^^ and Ramond'^*^^, Eq is uniquely se- 
lected among many GUT groups if we require i) every generation of quarks/lepton 
fields belongs to a single irreducible representation of the group, ii) the theory is 
automatically anomaly free, and iii) all the (phenomenological) Higgs fields neces- 
sitated for causing the symmetry breakings down to SU {3)c x C/(l)em fall in the 
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representations which can be suppUed by the fermion biUnears. Therefore it is very 
important to investigate GUT's from the viewpoint whether they are compatible 
or not with dynamical symmetry breaking. 

In this paper we study dynamical symmetry breaking in an Eq GUT model. 
The reason why we adopt Eq is its unique property stated above. In particular 
the third point implies the possibility that all the Higgs fields necessary for the 
symmetry breakings can be formed dynamically as fermion bound states. The 
unified gauge coupling constant of Eq suggested by the present experimental data, 
however, seems not large enough to break the Eq symmetry itself, and so we expect 
that some strong gauge interaction yet other than the Eq one exists and gives a 
primary driving force for the Eq symmetry breaking. But we still have no definite 
idea about that gauge interaction beyond Eq. So we assume in this paper that the 
strong gauge interaction is effectively treated as a Nambu-Jona-Lasinio type 4- 
fermi inter action [4]. We include all possible E'e-invariant 4-fermi interaction terms 
that can contribute to the formation of scalar Higgs fields. We, however, restrict 
ourselves to the model of a single generation of quarks/leptons, with the hope 
that the Higgs fields are all supplied as bound states of mainly a single generation 
of fermions. Following the usual procedure we introduce Higgs fields as auxiliary 
fields. We analyze the effective potential to find the the patterns of dynamical 
symmetry breaking realized in this model, and see whether the desirable symmetry 
breaking patterns emerge or not. 

This paper is organized as follows. In Sect. 2 we present the model which we 
study in this paper and give the effective potential of the auxiliary Higgs fields. 
The analysis of the effective potential is performed analytically in Sect. 3. In a 
special case in which the 4-fermi interaction is present only in Eq 27 channel, 
a complete analysis is possible and is given there. Otherwise, however, such a 
direct analysis becomes almost impossible and we perform a simplified analysis 
assuming that Michel's conjecture concerning the potential minimum holds. The 
symmetry breaking patterns found this way are actually confirmed to be correct 
by numerical analysis performed in Sect.4. Sect. 5 is devoted to the summary and 
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conclusion. Three appendices are supplemented; by using the spinor representation 
of 5*0(10) presented in Appendix A, representations 27 and 351 of Eq and an 
invariant tensor with three 27 indices are explicitly constructed in Appendix B; 
definition of maximal little groups which appears in Michel's conjecture is presented 
in Appendix C. 



2. The Model 



As explained in the Introduction, we consider Nambu-Jona-Lasino type model 
with a single generation of left-handed fermions, ip = (iPa) = 1, ■ ■ ■ , 27), be- 
longing to 27 representation of Eq. The Langrangian is given in the most generic 
form as follows: 



£, — £,Q + £int 

Co = - igA^)^^ + tr{Fi"'F^,) (2.1) 

Ant = + G'27|V'^C^|'27 + G351sl^^C^|'351s + GgsU I^^C^l'sSU ■ 



In this expression C denotes the charge conjugation matrix of Lorentz spinor so that 
ip'^Ct/j is Lorentz scalar, Gj's denote couphng constants, subscripts such as 27 mean 
the projection into the denoted irreducible component of Eq constructed with the 
fermion bilinear: 27x27 = 27+351s+351a- The absolute squares are understood 
to denote E'e-invariant contractions between those irreducible components and their 
complex conjugates. For the present case of single generation fermions the fermion 
bilinear ip^CipB is symmetric with respect to the indices A and B, and so the 
last anti-symmetric component (i/^^^Ci/jb) 351 a vanishes identically. Henceforth 351 
without subscript always denotes 351s. 

Now we introduce auxiliary fields (-^^27/351)^-6 standing for —{'ip'^Cip 3)27/ 351 
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and rewrite the interaction part i2int into 



Ant = -{(^IC^i?)27(i^27)'^'' + h.c} - Mlil{Hl^H27) 

-{(V^5CV'B)35i(//35i)''^ + h.c} - M|5i tr(i/|5ii/35l) , (2.2) 



27/351 — r< 

^27/351 

We evaluate only the fermion one-loop diagram for our effective potential. That is 
formally the leading term in 1/Ng expansion if we introduce Ng copies of our single 
generation of fermions. We neglect the Eq gauge interaction since it is expected to 
be weak. Then the 1-loop effective potential of H is given by 

(t>{H)^4>Q{H)+(t>i{H) 

d^. . . (2-3) 



MH)^-a[ ,l^lndet(MtM-p2) 
M = 2{H27 + H351) = 2H . 



Here J^d'^p denotes that the integral over p is defined with an ultraviolet cutoff 
\pe\ ^ ^ after making the Wick rotation to Euclidean momentum p — > pg. If 
the coupling constants are large enough this potential has a minimum away from 
the symmetric point H — and the £^6 symmetry is dynamically broken. We can 
determine the direction of the symmetry breaking by searching a minimum of this 
potential. 
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3. Analysis of the Effective Potential 



3.1. Case of 27 interaction only 

We first consider the simplest case in which only the 27 part of the 4-fermi 
interaction is present; namely, 

G27 ^ 0, 6-351 = . (3.1) 

Then clearly the Higgs vacuum expectation value (VEV) can appear only in the 
27 component: 

Hz^i = 0, H = H27 . (3.2) 

As explained in Appendix B, 27 representation of Eq is decomposed into 1 + 
16 + 10 under the 5*0(10) subgroup and the 27 Higgs field H — H27 is expressed 
in a 'vector' notation as 



V 



Ha 

\Hm) 



(a = 1,2,---,16) (3.3) 
(M=l,2,---,10) 



where the subscripts 0, a and M stand for the 5*0(10) singlet, 16 spinor and 
10 vector representations, respectively. V is embedded into the 27 x 27 matrix 
H — H27 by the help of the invariant tensor P"^^*-^ carrying three 27 indices, 
which is explicitly given in Appendix B: 

jjAB ^ ^ABCy^ (3 4) 



(See Eq.(B.ll) for the explicit form of this matrix.) 
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Any 27 can be Eq rotated into the following 'standard' reduced-form: 



\ / 



R + il 



m 







} SO (10) singlet (real) 
}SO{10) 16 

}the first component of 5'O(10) 10 (complex) 
}tlie second component of 5'O(10) 10 (real) 
}the third to tenth components of 5*0 (10) 10 



(3.5) 



{vq, R, 1 ,171 are all real) . 



This is seen as follows. First, starting from a generic form (3.3) of V, the spinor 
component can be rotated away by using the Eq rotation freedom with the 
spinor parameter e in (B.2). Next wc note that the vector component Hm actually 
stands for two 5*0(10) irreducible 10 vectors, the real and imaginary parts. So, 
using the S'O(IO) rotation freedom, we can make one of the two 10 vectors, say 
the imaginary parts, to have only the first component, and finally, by using the 
remaining SO (9) rotation freedom, we can make the real part 10 to have only 
the first two components. Alternatively, one can also convince the validity of the 
statement as follows: the problem is whether any 27 V can be written in the 
form gVo with g E Eq by using the standard form Vq in (3.5). Note as for g in 
this expression that only the right quotient Eq/S0{8) part is effective since Vq is 
invariant under 5*0(8). So the g part is parameterized by 78 — 28 = 50 parameters, 
and hence gVo spans a 50+4 = 54 dimensional space. But it is the same dimensions 
as the whole 27 complex vector V does. 

This standard form (3.5) has four parameters and implies that there exist four 
E'e-invariants which can be constructed by 27 representation V alone. They can 
easily be found and are given as follows: 



X = vW 



(3.6) 
(3.7) 



VaVeVc ^-iid its complex conjugate Y* 
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(3.8) 



These invariants X, Y, Y* and Z are expressed in terms of the four parameters in 
(3.5) as 

X = vl + R^ + rr? + I^ 

Y = vo{R^ + rr? - + 2iRI), Y* = vq{R^ + rr? - P - 2iRI) (3.9) 
Z = i^{R^ + rr? + Pf - APrr?^ + 4.vl{R^ + rr? + P) . 

Since the 1-loop effective potential (t){H) is invariant under Eq, (p{H) can be 
expressed in terms of the invariants, X, Y {Y*) and Z alone. The effective potential 
(2.3) now reads 



(f){H = M/2) = MI7- tr(MtM) - 4 / -f^ lndet(MtM - p^) 
4 J i{27ip 



1=1 ^ ' i=\ 



(3.10) 



where A^'s are real positive eigenvalues of M't'M, given by the roots of the following 
equations: 

Af - 8XA2 + leX^A, - 64FF* =0 for i = 1, 2, 3 (3.11) 
Af - 4X A| + 4ZAi - 16FF* = for i = 4, • • • , 27. (3. 12) 

These equations depend on only three quantities of the four invariants: 
X = -yg + + + 72 

yy* = t;2{(i?2+^2^j2)2_4j2^2| ^^^^^^ 
Z = |(i?2 + ^2^/2)2 _4j2^2| ^4^2^^2^^2^ j2^ _ 



For convenience, we re-parameterize these three quantities as follows: 



a = VQ X — a + b 

b = R^^rr?^P ^ yy* = a(&2 - 4c2) (3.14) 
c = Im Z^{b^ - Ac^) + Aah 

The three roots of Eq.(3.11) cannot explicitly be expressed in terms of a, h and c, 
but the twenty-four roots of Eq.(3.12) are given by 

Aj = 4a (i = 4...11) 

Ai = 26 + 4c (i = 12...19) (3.15) 
Ai = 26 - 4c (i = 20 . . . 27) . 

Then, inserting these exphcit expressions for the roots Aj (i = 4, • • • , 27) and using 
an identity 

(Ai + y){\2 + y){\i + y)=y'' + SXy'' + IGX^y + 64YY* (3.16) 

following from Eq.(3.11) for the implicit roots Ai^2,3, the effective potential (3.10) 
reduces to 

cP{H) = cP{a,b,c) 

= 10M|7(a + b)-^j ydy [\n{y^ + SXy'' + IGX'^y + (3.17) 
+ 8 ln(4a + y) + 8 \n{2b + 4c + y) + 8 ln(26 -Ac + y) 

We now look for the stationary point of the effective potential (p{H) — (f){a, b, c). 
Taking account that X is independent of c and d{YY*) / dc — —Sac, the derivative 
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of the potential with respect to the parameter c is given by 



^ [-64-8ac/(A2) + 8-4{c/(26 + 4c) -^(26-4c)}] (3.18) 



dc An 
with functions / and g defined by 



Jo 



y 



g{x) ^ dy 
Jo 



+ 8Xy2 + lex'^y + 6AYY* 



— A — x in 



(3.19) 



y + X X 
Note that f{x) is positive for x > since /(O) = and 



•^'^^^ " x3 + 8Xx2 + i6X2x + 64yy* ^ ° ^^"^^^ 



because X > and yy* > by definition (3.14). On the other hand, g{x) is seen 
to be a monotonously decreasing function of x since 

g'{x) = - r < . (3.21) 

Jo {y + xf ^ ' 

Taking account also that c is bounded (|c| < \h\/2) by definition (3.14), we find 

that 

sgn(|^) = sgn(c) . (3.22) 
This shows that 0(a, 6, c) has an absolute minimum at 

c = (3.23) 

in the defining region |c| < \h\/2. 
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Next consider the derivatives of 0(a, b, c) with respect to a and 6; at c = they 
are given respectively by 



90 
da 



IOMI7 - 



47r2 



ydy 



8|/2 + 32Xy 



+ 64-b^ f{A^) + S-4g{4a) 



(3.24) 



86 n 1 

db ^' 47r2 



A2 



8|/2 + 32X|/ 



+ 64-2abf{A^) + 8-4g{2b) 



(3.25) 



Stationarity requirement d(f)/da = and d(f)/db = 0, or only the difference d(j)/da — 
d4>/db = 0, leads to the following condition: 



26(2a-6)/(A2) = ^(4a)-^(26) . 



(3.26) 



Since a, 6 > by definition and /(A^) > as mentioned above, the sign of the LHS 
is sgn(2a — b) while the sign of the RHS is opposite, — sgn(2a — 6), since g{x) is a 
monotonously decreasing function. Therefore 0(a, b, c) can have a minimum only 
at 



2a 



(3.27) 



Under the conditions (3.23) and (3.27), the 27 vector (3.5) now takes the form: 



V 



X } so (10) singlet (real) 





V2v 



}SO{W) 16 

}the first component of 5*0(10) 10 (real) 
V ^ / }the second to tenth components of 50(10) 10 



(3.28) 



For the VEV of this form, the 27 eigenvalues Aj of the fermion squared mass 
matrix MtM, determined by (3.11) and (3.12), now become explicit and show an 
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interesting degeneracy: one 16v^ and twenty six Av'^^s. Then the potential (3.10) 
is given by 

1 r 1 

(l){H) = 3QMijV^ [ln(16v2 + y) + 26 ln{4v'^ + 7/)J , (3.29) 



the stationarity of which determines the magnitude of the VEV v: 

1 

47^2 



3OMI7 - (16 g{16v'^) + 26-4 g{Av^)) . (3.30) 



The critical coupling Gf^ — l/(-^27)^ beyond which non-zero VEV is reahzed is 
found by taking — > 0: 

Gf7 = ^ . (3.31) 

Note that on this vacuum one fermion has mass 4v and all the other 26 fermions 
have a degenerate mass 2v. This implies that the original fermion multiplet 27 
splits into 1 + 26. This branching pattern indicates that the symmetry breaking 
realized in this case is 

Ed — > F4 . (3.32) 

This can also be confirmed by calculating branching pattern of the Eq gauge boson 
masses on this vacuum; 78 26 (massive) + 52(massless), where 52 is massless 
gauge bosons of the unbroken F4 group. 

3.2. Case of 351 interaction only 

Next we consider the case in which only the 351 part of the 4-fermi interaction 
is present; namely, 

G351 7^ 0, ^27 = , (3.33) 
in which case the Higgs VEV appear only in the 351 component: 

H27 = 0, Hs5i . (3.34) 

Contrary to the previous case, there are many Eq invariants and it is almost im- 
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possible to perform the same kind of analysis as in the previous subsection. 

Now we invoke Michel's conjecture[l,ll], which claims that the following state- 
ment holds when the potential system contains only a real irreducible representa- 
tion of scalar fields, or a self-conjugate pair of a complex irreducible representations: 

The group symmetry can break down only to one of the maximal little groups of 
the representation considered. (See Appendix C for the definition of maximal little 
groups.) 

For illustration, let us apply this conjecture to the previous case; there, the 
fields appearing in the potential are (auxiliary) Higgs fields of 27 and its conjugate 
27. Then the maximal little groups arc 5*0(10) and F4. We in fact found the 
symmetry breaking Eq —>■ F/^ in the above, so that the conjecture was actually 
correct. 

In the present case of 351 Higgs fields, the maximal little groups are 5'O(10), 
F4, SpiS), G2, SU{3) and SU{2)®SU{A). We assume that the conjecture holds in 
this case also. Then what we have to do is to calculate the effective potential for 
each possibility of the Higgs VEV's in the maximal little group directions and to 
compare the minimum values of the potentials to see which possibility is actually 
realized. 



He Eg 


decomposition under H 


Fa 


27= 


=1+26 


Sp{8) 


27= 


=27 


G2 


27= 


=27 


SU{3) 


27= 


=27 


SU{2) ® SU{4) 


27= 


= (2,6) + (l, 15) 



Table 1. Decomposition of Eq 27 under the maximal groups H. 
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Let us first determine the form of the VEV for each case of the maximal httle 
groups. From Table 1, we see that the fermion 27 of Eq is again 27 under 5*^(8), 
G2 and SU (3) but the latter is a real representation. Taking also account that Eq 
351 is constructed from 270 27 symmetrically, we see that the trace part of i^asi 
is a singlet under those groups. Therefore, in a suitable basis, the VEV takes the 
following form for the cases of those groups: 



H351 = v(®l27 {Sp{8), G2, SU{3)) , (3.35) 

where 1„ denotes n x n identity matrix. 

Similarly, in case of5'C/(2)(8)5'C/(4), a singlet under this group which we get 
from a symmetric product of 27 x 27 comes from the component (1, 15) x (1, 15) 
and hence, in a basis, 

/^^®ll5 A 

Hs5i = (SU(2) SU(A)) . (3.36) 

V ® I12 / 



It is a bit more complicated to get the form of VEV in cases of SO (10) and 
F4, since, as is seen from Table 1, we get two singlets from the symmetric product 



27x27for eachcaseof 5'O(10) and F4. In view of 27(8)27|sym = 27+351, we see 
that one singlet is in H27 and the other is in i^asi. We can find the form of 5'O(10) 
and F4 singlets in 27; from (3.4), (3.5) and (3.28), they are given, respectively, by 



(g) li7 

^27=1 , 1 (SOm). (3.37) 

V (g) lio 




H27^{ (^4). (3.38) 

We see from (3.37) that the SO{10) singlet in H27 comes solely from 10 x 10. 
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Therefore the 5*0(10) singlet in ifasi must be the other singlet made from 1x1: 



^351= r I (50(10)). (3.39) 

0® l26 

(This can also be seen directly from (B.IO) in Appendix B.) In case of F4, the 
two F4 singlets come from 26 x 26 and 1x1, one combination of which is (3.38) 
contained in H27. Since 351 is orthogonal to 27, tr H^jHs^i = 0, the F4 singlet 
in i?35i should thus have the form: 

. 13^; 

H35i^{ I {F4) . (3.40) 



V (g) 1 



26 



Now that we have found the form of VEV's, we can calculate the minimum 

value of the potential by substituting those matrix into (2.3) for each case and 
compare their minimum values to find the direction of symmetry breaking. First 
we define a function 

F{v) ^ - 4 ln(4^2 - p') . (3.41) 

Then the potential value in the directions 5p(8), G2 and SU{3) is commonly given 

by 

<Psp = 27F{v) . (3.42) 
In the directions SU{2) 5C/(4), 50(10) and F4, it is given respectively by 

cf>ss = l5Fiv) , (3.43) 
0SO = F{v) , (3.44) 
(f)F4 = 26F{v) + F{13v) . (3.45) 

When symmetry breaking occurs, F{v) has a minimum at a certain point vq and 
takes a negative value there. The potentials (3.42), (3.43) and (3.44) take their 
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minima at the same point vq and hence we immediately see for the minimal values 

(l>Sp < (l>SS < (t>SO ■ (3.46) 

The minimum of (3.45) is realized at a certain point vi, which is different from the 
minimum point vq of F{v), so that 

F{vq) < F{vi) and F{vq) < F(13vi) , (3.47) 

and hence 

(t>Sp < (t>FA . (3.48) 
We thus find that the symmetry breaking in this pure 351 interaction case is 

Eq — > Sp{8) or G2 or SU{3) . (3.49) 

These three group cases cannot be distinguished in the present approximation in 
which only the fermion one-loop vacuum energy is counted, since the fermions get 
quite the same masses for those three breakings. This degeneracy will be hfted if 
the vacuum energy due to gauge boson loops are taken into account. 

3.3. GENERAL CASE 

Finally in this section we study the general case in which both 27 and 351 4- 
fermi interactions are present. Strictly speaking, Michel's conjecture is inapplicable 
to this general case since there appear two fields of different representations 27 and 
351 in the potential. Nevertheless we assume that this conjecture still holds and 
determine the symmetry breaking pattern in this case also using the same analysis 
method as in the previous subsection. 
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Candidate groups arc the same as the 351 case: 5*0(10), -F4, 5*^(8), G2, SU{3) 
and SU{2) ® SU{A). Of these groups Sp{8), G2, SU{3) and SU{2) ® SU{4) have 
their singlet only in 351 of Eq, and so the VEV's and potentials are the same as 
in the previous subsection. Therefore 

<Psp < <Pss (3.50) 

is always realized. 

On the other hand 5*0(10) and F4 have their singlets in both 27 and 351 
representations of Eq. Their singlets in 351 arc contained in the form (3.39) and 
(3.40) and those in in 27 are in the form (3.37) and (3.38), for the 5*0(10) and 
F4 cases, respectively. Thus the general forms of VEV's for these group cases are 
given respectively by 



H 



(V \ 

0® I16 



(-50(10)) (3.51) 



, 13^; - 2V 

H=\ (F4) . (3.52) 

By using (3.51), the potential corresponding to 5*0(10) breaking (= 0^q) is 

= (^It - ^Isi)^' + F(y) + lOi^(^) (3.53) 
and by using (3.52), that of F4 breaking (= 0^4) is 

0'^4 = 30(M|7 - M|5i)v2 + F{Uv - 2V) + 2QF{v + V) . (3.54) 

Now let us compare (pspj4>'so their minimum points. First of all we 

clearly see the relation 

(pSp = 0F4 < (f>'so when MI7 = MI51 . (3.55) 

In view of this, we study the potential in two cases of (a) M27 > M^^-^^ and (b) 
M27 < M^r,i, separately. 
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(a) MI7 > M3251 

In this case we have from (3.53) 

'P'so = (Mir - Ml^)V^ + F{V) + 10F(^;) 

> llF(^;o) (3.56) 

> 21F{vq) = (t>sp , 

and from (3.54) 

(f)'p,^ = 30(M|7 - M|5i)\/2 ^ p^i^^ _ 2V) + 26F{v + V) 

> F{13v - 2V) + 2QF{v + V) (3.57) 

> 21F{vq) = <t>sp . 

Hence we conclude that the symmetry breaking pattern in this case is given by 

Eq — > Sp{8) or G2 or SU{3) (3.58) 

(b) MI7 < Mi,, 

Taking account that 4>sp does not depend on M27, we first study the derivative 
of (p'p^ with respect to MI^, with M^,-, kept fixed. The arguments V and v of 0'^4 
arc set equal to the values realizing the stationary point of (p'p^ and so they depend 
on MI7. 

5(^'^4 (stationary point) _ 3 , 9V d(p'p^ dv d(t)'p^ 

dM^j ~ dM^j dV 9M|7 dv (3.59) 

= 30V^ > . 

This implies that the minimum value of (f)'p^ is monotonously increasing as a func- 
tion of M27, and hence together with (3.55) that 

<P'f4 < (t>Sp (3.60) 

in this region Mfy < Mgg^. 
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Next we compare (p'p^ and (p'gQ- In the limiting region 

MI51 > MI7 namely G351 < G27 ^ 00 , (3.61) 

the system is the same as that where there is only 27 4-fermi interaction and there, 
as we know, the F4 vacuum is the lowest one: 

4'W < 4''so ■ (3-62) 

On the other hand the relation (3.55) implies that the same inequality holds even 
in the region M27 ~ ^351- This strongly suggests that the inequality (3.62) holds 
for the whole region M27 < -/Wgsi- We assume this holds. Then, together with 
(3.60), we find that the symmetry breaking pattern in this coupling region is 

Eq — ^ F4 . (3.63) 

The discussion in this subsection is very incomplete by two reasons. Firstly, this 
general case is outside the scope of Michel's conjecture. Secondly, the Eq.(3.62) was 
not proved for the whole region of M27 < M^^-^. Nevertheless it suggests a simple 
symmetry breaking pattern; it is either Eq F4 or Eq Sp{8) or G2 or SU (3) 
depending on whether the 27 interaction G27 is larger or smaller than the 351 
interaction (^351, respectively. 
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4. Numerical Analysis 



In order to confirm the symmetry breaking pattern suggested by the analysis 
in the previous section, we numerically search the minimum of the potential (2.3) 
and calculate a fermion mass spectrum and gauge boson mass spectrum at that 
point. 

4.1. Algorithm 

We present in this subsection an algorithm for searching the stationary point 
Hst- = 0. The idea is essentially to apply the Newton method to the 

derivative d(f){H)/dH'^ since we want a zero point of this function. 

First of all we note: 

1) 4>{H) is a function of 378 x 2 variables as is a 27 x 27 symmetric and 
complex matrix. 

2) d(l){H)/dH^ = V{H) is a gradient of in the 756 dimension space, which 
can be written down in a closed matrix form: 



3) On the contrary we have no such a simple analytic expression for the second 
derivative of 4>{H). 

We now outline how the iteration method goes for searching the minimum. (We 
assume in the following for simplicity that (f){H) is concave in the considered re- 
gion.) 

i) We take randomly a starting H = Hq, and calculate the gradient V{Hq) — 



V{H) 




(M|7 - Ml^)H27 + Ml^H 



(4.1) 




offiHo) there. 
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ii) To find the next point wliicli is nearer to tlie stationary point, we consider the 
potential function in a cross section in the gradient direction; namely, 

we consider the following function of one real parameter t: 

f{t) = <P{HQ + Vt) . (4.2) 

If we find a zero of the first derivative function 



. N df(t) 



(4.3) 



at t = to, then H = Hq + Vto will be the lowest point of 4>{H) in this cross 
section. 

iii) Starting from t — 0, the Newton method applied to this function g{t) gives 
at the first iteration step 

U^-^ (4 4) 

as a nearer point to the zero to of g{t). We do not continue this Newton's 
iteration any further since even if to is found more exactly the point H — 
Hq + Vto is merely the lowest point of (f){H) inside this cross section. So we 
adopt Hi — Hq + Vti as a nearer point to a true stationary point of (f){H). 

iv) If ti is already small enough we consider Hq is a stationary point. Otherwise 
we take Hi = Hq + Vti as Hq in the step i) and repeat the procedure. 

What we get by this iteration procedure is, logically speaking, not a minimum 
point but a stationary point. But, in practice in this calculation, we actually 
obtained a minimum although it may not be a global minimum. 
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4.2. Result 

We have run the above procedure for searching the stationary point for the 
potential with various sets of parameters, M'^^-^/M2'j and A^/M27; more cxpUcitly, 
we have swept the region < Mfg^/Mf^ < 10^ and 40 < AV^27 ^ 10^- (Note 
that A^/M27 = TT^ is the critical value for the symmetry breaking in the pure 27 
interaction case.) We have stored in total about 10^ data of the stationary points 
Hst for the potentials with the parameters in this region, in particular, in the region 
10-3 < Mfg^/Ml^ < 300 and AV^It = 40, 100 in detail. 

Using the obtained stationary point data i?st, we have calculated fermion 
masses and gauge boson masses on those vacua. Fermion masses are calculated as 
eigenvalues of the squared mass matrix H^^Hgi and those of gauge bosons are as 
eigenvalues of the squared mass matrix G = {Gab) = (ti'[Ta-?/j^+-f/^]tT^])(tr[T^*i^st+ 
HstT^]) where TaS are Eq generators in the 27 representation, whose explicit form 
is given in Appendix B. We can judge the symmetry breaking pattern from those 
mass spectra for each case. 

The result of our numerical calculation is summarized as follows. 

1) When M27 < M351, namely, 27 4-fermi interaction is dominant, we found in 
every case of our search the following. 26 fermions has a degenerate mass and 
the rest one fermion has another mass. On the other hand, 52 gauge bosons 
are massless and the rest 26 has a degenerate non-zero mass. All these clearly 
imply that the symmetry breaking pattern in this coupling region is 

Eq F4. (4.5) 

This completely agrees with the result obtained in the previous section, de- 
spite that the latter was based on a bit non-rigorous arguments. 

2) When M27 > M351, namely, 351 4-fermi interaction is dominant, we found 
in every case the following. All of the 27 fermions has a degenerate mass 
while the gauge bosons become all massive but not degenerate at all. The 
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degenerate fermion spectrum implies that the symmetry breaking pattern in 
this case is 

Eq — > Sp(8) or G2 or SU(3) , (4.6) 

agreeing again with the result of the previous analysis. But it seems strange 
why all the gauge bosons are massive and non-degenerate. If there remains 
some symmetry, the corresponding gauge bosons should remain massless and 
the spectrum should show some multiplet structure. The reason why this 
strange thing happens is in the particular nature in this breaking; namely, in 
this case, the three different vacua with symmetries Sp{S), G2 and SU (3) are 
degenerate. They place at different points in the potential but realizes the 
same stationary value. Then, if there is a path connecting these three points 
through which the potential is fiat (or almost fiat within the calculation 
error), all the points on the path realize the same stationary values but has 
no symmetries at all. Nevertheless, the fermion mass degeneracy is still 
realized since the present effective potential counts only the fermion vacuum 
energy and the degeneracy of the potential value along the path means the 
fermion mass degeneracy. All the stationary points we found are such points 
on the path. This is our interpretation, but we confirmed this by examining 
the potential values realized by our stationary points. They all coincided 
with (psp — 27/(i'o) which we obtained analytically in the previous section 
by using Michel's conjecture. 





M27 < M351 


M27 > M351 


fermion mass 


1+26 


27 


gauge bosson mass 


52 (massless) + 26 (massive) 


1 (massive) X 78 



Table 2. Mass spectra found numerically for the cases M27 < M351 and 
M27 < M351. 
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5. Summary and Conclusion 



We have analyzed an Eq GUT model of a single generation of fermions with 
strong 4-fermi interactions. The Eq symmetry is found to be broken spontaneously 
either to F4 or to 5*^(8) or G2 or SU (3) depending on which of the 4-fermi coupling 
constants G27 and G351 in the 27/351 channels is stronger than the other. 

In these symmetry breakings, the fermions turn to belong to real representa- 
tions of the residual symmetry and all of them get acquire non-vanishing masses. 
Since these masses arc necessarily of the order of the GUT symmetry breaking 
(~ lO^^^^^GeV), the present model as it stands, unfortunately, turns out to be un- 
realistic as a GUT model. The quarks and leptons belong to a chiral representation 
of the standard gauge group and should remain massless at the GUT scale. 

We can easily understand the reason why all the fermions get non-vanishing 
masses in the present model. As mentioned before, our effective potential counts 
only the fermion one-loop vacuum energy. But the fermion vacuum energy es- 
sentially comes from the energy of Dirac's negative energy sea and hence is neg- 
ative. So, the more massive the fermions become, the more the vacuum energy 
is lowered. Therefore the desirable symmetry breaking patterns, such as down to 
SU{3) X SU{2) X C/(l) under which the fermions are chiral and remain massless, 
are necessarily disfavorable energetically. 

This indicates that the vacuum energy coming from bosons should play a cen- 
tral role in order for the present model to produce desirable symmetry break- 
ing patterns. Indeed, Harvey[12] once considered the Eq symmetry breaking in a 
Goleman- Weinberg like spontaneous symmetry breaking scenario and found that 
Eq is broken down to SO (10). There the main part of the potential in fact came 
from the gauge boson loop contribution. 

Alternatively, there may be another possibility if we change the fermion content 
of the model. For instance[13] , we can regard the three generations of quarks/leptons| 
as merely survivals from GUT world where n+3 generations and n anti-generations 
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of fermions exist. Then, when the dynamical GUT symmetry breaking occurs, a 
variety of mixing can generally occur among those fermions, and n generations 
of fermions as a net number can get acquire O(Mgut) masses leaving the usual 
quarks and leptons massless. Since there are fermions which aquire the masses in 
this case, there is a possiblity that small contributions of the gauge boson loop 
may be sufficient to realize such desirable breaking down to chiral type symmetry. 
This type of scenario is very interesting also from the viewpoint of the origin of 
Cabibbo-Kobayashi-Maskawa mixing as well as of the stabihty of proton. 
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APPENDIX A. 50(10) 7-Matrices 

For any SO{2n), the 7-matrices ^^Fm (M = 1, 2, • • • , 2n) satisfying ^Tiv+I 
2np^ 2np^ ^ 2(5Miv take the form 

where and are 2 a' "^-component Wcyl spinors with chiralities '^'^T2n+l = +1 
and —1, respectively. The superscript on the left shoulder indicates the dimension 
2n of SO{2n). cr's are the 7-matrices in the Weyl spinor basis. 
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Totally anti-symmetric multi-indexed 7-matrices '^^TMiM2-Mk ^-^e defined by 



2n 



1 



^M.M^-Mu = ^ ( • ' • ^"Tm, + (anti-symmetrization)) 



( / 2n. 










2rb— 

<^MiM2--Mk 

(7MiM2--Mk 





for k — even 



for k — odd 



2n 



^ ( ^'Vmi ^"(jif, ^Vms ^"VJ^^ • • • ^Vj^ + (anti-symmetrization)' 

^ ( ^"-^Mi ^"-^Ma ^Vj^^3 ^'Vm4 • • • ^'V^J, + (anti-symmetrization)) 

(A.2) 

The S0{2n) generators Tm at satisfying [T^iv, ^r-l] — —'i{^NKTML+ (anti-symmetrization) )| 



are expressed in this spinor representation by the matrix 

2rv - 1 2np ^ ( '^'^MN 

^MN = 777 i MN = 9^ 

2^ 2z \^ 2r^MiV 

The charge conjugation matrix ^"C exists such that 

2«^^2np^2r^-l ^^2npT 

2r^t 2r^ ^ 1 ^ 2r^T ^ ^ 2r^ ^.^j^ ^ = COS ^ + sin ^ 



(A.3) 



(A.4) 



for either choice of 77 — ±1, where the superscript T denotes transposed. Hence- 
forth we always choose 77 = -1-1 for convenience. 

A.l. 50(6) 

We first construct 5'(9(6) 7-matrices in the Weyl spinor basis: it is convenient 
to take the 4x4 matrices as 

(y-m = [ (^1=1,2,3, cri+3=4,5,6 ) 

{^^i)aP^^i^ocP + 5% (A.5) 



where SaP'yd is rank-4 totally anti-symmetric tensor and (J^s is multi-index anti- 
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symmetric Kroncckcr's delta defined by = ^^^^j — ^a^^- The index i here, 
running over 1, 2, 3, will correspond to the color index of SU{2>) C SO{Q). 

These (m = 1, 2, • • • , 6) possess the following properties: 

^c"m = ~ ^CTm ( ant i- symmetric) 
1 

( ^'^^) a(3 ^ ~ 2 ^"'^^'^ ^ ^'^'"'^ ^ (anti-selfduality) 

An 5'0(6)-vector Vm is eqiuvalent to a rank-2 antisymmetric tensor V[q,^] of SU{A)] 
they are related with each other via 

Ig ^ ( ^ Q \ \ '^^ 

V[af}] = ^( '7m)apym ^ Ki = -(^^ C^mj ^[ap] ■ (A.7) 

Decomposition of the 5*0(6) vector Vm into 3 + 3* under the color group SU{2)) C 
SU{A) ~ 50(6) is given by 

1 

3' 



3: V[i4] = ^(yi-zl^i+3) ■ 



(A.8) 

The 50(6) generators are given by the general expression (A. 3), which defines 
^(Jmn and ^mn- Then the 15 matrices ^(Jmn ("^5 = 1? " ' ' > 6) together with a unit 
matrix span a complete set of 4 x 4 matrices and satisfy the following completeness 
relation: 

^(l)J(l)/ + ^('^mn)J( Wn)/ ^ ^i^} . (A.9) 

The charge conjugation matrix ^ defined generally in (A.4) is now given by 
with Im denoting m x m unit matrix. 
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A.2. S0{4) 



SO (A) 7-matrices (fi = 7, 8, 9, 0) in the Weyl spinor basis are 2x2 matrices 
which we take as follows: 

^CTfx ^ {-icri,-ia2,-icr3, 12) (A. 11) 

with cTi,2,3 being the Pauh matrices. Then SO{A) generators ^^i^ defined by (A. 3) 
split into 3 + 3 generators of SU{2)j^ x SU{2)ji ~ SO (A): 



1/1 \ /o 



(A.12) 



The charge conjugation matrix 'kj is given by 



t7=('^' ° )=-t7T. (A.13) 

i(72 



A.3. 50(10) 

5*0(10) 7-matrices ^^r^vf (M = 1, 2, ■ ■ ■ , 9, 0) are constructed by a tensor prod- 
uct of the 50(6) and 50(4) 7-matrices as follows: 

1. ( ^T„,0^T5 ior M^m^ 1,2,- ■■,6 . ( U 

= { , , with T5 = 

( 14(8) % for M = = 7, 8, 9, \ -I2 

(A.14) 

Then the 7-matrices in the Weyl basis, ctm, for which we omit the superscript 10 
implying 50(10) for notational simphcity, are 8x8 matrices taking the following 
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form: 



0'M=m -fit-. ^ ' (^M=ii - I ^ ^ 4 + 

12 y \ I4 V], 

(A.15) 

The charge conjugation matrix takes the form 

C 

iOC^(^^^^ [ 1 (A.16) 



C 



where C is 8 x 8 matrix given by 



, -I4 ® i(72 , rr 1 + 

C = ( I = = = . (A.17) 

I4 ® ■^0-2 

Prom ^^'^ = it follows that the matrices ^'^m are symmetric, and 

(t) (t) — 

so are Ca)^ and o^C. Similarly we see that Ca MiM-z-'-Ms and crMiM2---MsC are 

symmetric. Since they are selfdual, 

CctMi-Mb = -^^i^Mi-M^Ni-NsC'^Ni-N^, (A. 18) 

they give loC's/S = 126 symmetric matrices, and hence, together with the ten 
Ca\^ matrices, span a complete set in the space of 16 x 16 symmetric matrices; 
the completeness relation reads 



2-4 



1 1 



2 

(A.19) 

because of the normalization condition 

2""' tr (Ca\^aNC) = Smn 

^ ' (A.20) 

2~'^tr {CaMx-Mz(yN^-N^,C) = (^Mi'.'.'.Ms + i^Mi-MsNi-Ns ■ 

In the same way 10C3 = 120 matrices CCTM1M2M3 (or, (TM1M2M3C) turn out to give 
a complete set of 16 x 16 anti-symmetric matrices so that 
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APPENDIX B. Some Representations of 



The Eq algebra is most easily expressed by refering to its maximal subgroup 
5*0 (10) X U{1). The generators are given by 16 SO (10) Weyl-spinor generators 
Ea (« = Ij ■ ■ ■ J 16) and their complex conugates = (-Ea)^ in addition to the 45 



5*0(10) generators Tmn and one U{1) generator T. The algebra is given by 
= -i {SnkTml + SmlTnk - SmkTnl - ^nlTmk) , 



[15] 






(B.l) 



[Ea, E^] = --{crMN)aTMN + -^^^T . 



The simplest representation of Eq is 27 which is decomposed into I4 + I61 + 
10-2 under the maximal subgroup 5'O(10) x C/(l). (The suffices denote the value of 
U{1) charge 2-\/3T.) So the 27 representation can be denoted as = (■0o, '4>a, '4'm) 
with a and M being 50(10) (Weyl-)spinor and vector indices, respectively. The 



Eq generators act on this representation as 



[15] 



{OT^-OklTkl + e^E^ + 'e' e^) 



V 



e/3 







\OKL{(^KL)a + 



2^3 



UN — -^^^M 











J 


Uiv/ 



(B.2) 



To check that this representation for the Eq generators really satisfies the algebra 
(B.l), we need the following identities for the 50(10) 7-matrices: 

T T 
CumnC — — (J^s , CaMiMiMaM^C — (7M1M2M3M4, > 

(^MNCTR - CTk'^MN ^ i{5MK<7N - SnRCTm) , (B.3) 
l{Calf^aMC)sa - SIS^ = - ^(<7Miv)/(<TMiv)," . 
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The last identity follows from the Fierz transformation of the LHS: 



4! 



-j{C2iaMNC)^oi{(^K'^i(TMN(7^K)5 + -^^{C (7 Mi---MiC!)^g^{a kCT Mi---Ma^^k) 5 

(B.4) 



Tensor product of two 27 representations gives 



27 X 27 = 27s + 351s + 351a • (B.5) 
This implies that there is an invariant tensor F"^^^ which gives 27 from 27 x 27: 

= r^^^VfiV'c ■ (B.6) 

This r^^^ is found to be given by 

r ro^^ = 5mn 

r^^^ : totally symmetric I T^^"^ = ^(CcrJ^^)"'^ (B.7) 

L otherwise , 
j\ 

or equivalently, in terms of the components of ^ , 

= ^V^^Cal^^ + 2V'oV'M (B.8) 
= V2V^M(<::^(TLfV')" ■ 
To check that this ^ transforms correctly as 27, we need an identity: 

(e^Cal^lj)(^lj^Calrj) = -^{e^CaMV){^''Cal^^;) , (B.9) 

which follows from Fierzing cr^^ijj and rj and using cr\^aKcr\j — —8a\^ and a\^aK:i--Ks<^ 
0. 



The 351 can be represented by a symmetric tensor $ with two 27 indices A 
and B: 

4,$^ (16) 



^/2■ 

$"/3(i26) 



V2- 



-^^$"^(144) 

$^^(54) y 



(B.IO) 



where the argument in each entry denotes the dimension under SOiVS). The 
previous 27 representation ■^"^ can also be imbedded into a symmetric matrix 
using the invariant symmetric tensor Y^^*^: 



^ABC 



/ \ 



(B.ll) 



Since symmetric tensor product of two 27 is either 27 or 351, the 351 matrix 
^"^^ can be characterized as a general symmetric matrix which contains no 27 
components of the form (B.ll): therefore, the component should be traceless, 
^MM _ g. gj-^ould contain no 50(10) vector components, {crMC)ai3^'^^ = 
0; ^^/^ should be 7-tracclcss, ((TMC)a/?<^'^^ = 0. But, matter of course, 
these conditions are nothing but the requirements that those entries be irreducible 
representations under 5*0(10) as indicated in the arguments in (B.IO). 

APPENDIX C. Maximal Little Group 



A little group of a representation vector of a group G is defined by 



(C.l) 



This little group depends not only on the representation but also on the vector 
itself. 

Consider a single irreducible representation i? or a self-conjugate pair R + R* 
of a complex irredicible representation R. For this representation R, many little 
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groups appear as the vector (j) varies in the representation R with the length |0| (7^ 0) 
kept fixed. A httle group H is called maximal if there is no with little group 
satisfying G D D H. 

Some examples of Eq maximal little groups are given in the following Table. 



R 


Maximal httle groups 


78 


SU{Q) X [/(I), 50(10) X [/(I), SU(b) X SU{2) x [SU{2,)f x SU{2) x U{1) 


27 


50(10), F4 


351 


50(10), F4, 5p(8), O2, 5t/(3), 5t/(4) x SU{2) 



Table 3. Maximal httle groups for the representations it! (+-R*). 
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